arXiv:1501.06355v4 [astro-ph.CO] 23 Aug 2015 


Prepared for submission to JCAP 


Subdominant Dark Matter sterile 
neutrino resonant production in the 
light of Planck 


L. A. Popa, D. Tonoiu 

Institute of Space Science, 

Bucharest-Magurele, Ro-077125 Romania 

E-mail: lpopa@spacescience.ro, tonoiud@spacescience.ro 

Abstract. Few independent detections of a weak X-ray line at an energy of ~3.5 keV seen 
toward a number of astrophysical sites have been reported. If this signal will be confirmed 
to be the signature of decaying DM sterile neutrino with a mass of ~ 7.1 keV, then the 
cosmological observables should be consistent with its properties. 

In this paper we make a coupled treatment of the weak decoupling, primordial nucleosyn¬ 
thesis and photon decoupling epochs in the sterile neutrino resonant production scenario, 
including the extra radiation energy density via A^e// • We compute the radiation and matter 
perturbations including the full resonance sweep solution for Val^a flavor conversion 

in the expanding Universe. 

We show that the cosmological measurements are in agreement with subdominant Dark 
Matter sterile neutrino resonant production with following parameters (errors at 95% CL): 
mass =6.08 ± 3.22 keV, mixing angle sin^ 20 <5.61 lepton number per flavor 

Li = 1.23 ± 0.04 {Li = lO^Liy^) and sterile neutrino mass fraction < 0.078. 

Our results are in good agreement with the sterile neutrino resonant production parameters 
inferred in Ref. [62] from the linear large scale structure constraints to produce full Dark 
Matter density. 

Keywords: cosmology: cosmic microwave background, cosmological parameters, early uni¬ 
verse, neutrino 
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1 Introduction 

With about 25% of the total energy density of the Universe represented by Dark Matter 
(DM) and six basic cosmological parameters, the ACDM cosmological standard model is 
remarkably successful at reproducing the large-scale structure (LSS) of the Universe [1, 2]. 
Despite its successes on large-scales, this model produces too much power at small scales. 
In general, the observed structures are less abundant (“missing satellite problem”) [3-6], 
have cored density profiles (“core/cusp problem”) [7-9] and are too dense (“too-big-to-fail 
problem”) [10, 11] than those predicted by the ACDM model. Numerous attempts involving 
collisionless DM [4, 12-15] fail to resolve these problems, suggesting that DM may exhibit 
gravitational properties as well. 

The nature of DM remains however unknown and the explanation of its existence requires 
the revision of Standard Model (SM) of particle physics. In this regard, many extensions of 
SM have been considered [16-18] . 

Among the highly motivated candidates for DM is the right-handed sterile neutrino [19- 
22]. The minimal extension of SM by the addition of three right-handed sterile neutrinos (the 
i/MSM), each one corresponding to a SM active neutrino flavor, can simultaneously explain 
the active neutrino oscillations and the DM properties [23, 24]. Within this theory, the DM 
sterile neutrino masses are in the keV range and the mass of the lightest active neutrino is 
smaller than 10“^ eV. This excludes the degenerate mass spectra of the active neutrinos and 
hxes the absolute mass scale of the other two active neutrinos [25]. Moreover, the i/MSM 
can also explain the matter-antimatter asymmetry of the Universe [24]. 

In the original Dodelson-Widrow (DW) model [26], sterile neutrinos are produced in the 
early Universe through non-resonant oscillations (NRP) with active neutrinos in presence of 
negligible leptonic asymmetry. As the DW sterile neutrino is a Warm Dark Matter (WDM) 
particle with non-negligible velocity dispersion, this suppress the matter power spectrum 
bellow the free-streaming scale, that affects the DM distribution at small-scales. 

For NRP case, the universal Tremaine-Gunn bound [27] on sterile neutrino mass leads to 
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j^NRP ^ g ]^gY_ Lower limits on sterile neutrino mass have been placed from various 
observational probes. The combination of CMB measurements, LSS and Ly-a forest power 
spectra lead to >1.7 keV with a further improvement to > 3 keV when high- 

resolution Lyman-a power spectra are considered [28]. By modeling the SDSS Lyman-a 
forest flux spectrum, a lower limit of > 13 keV have been found [29, 30]. The same 

lower limit has been placed by requiring the number of subhalos in N-body simulation to be 
larger than the number of the observed dwarf spheroidal galaxies (dSphs) of the Milky Way 
[31]. Using the phase-space densities derived for dwarf satellite galaxies of the Milky Way a 
recent work [32] found a lower limit of DW sterile neutrino mass of >2.5 keV. 

Sterile neutrino with a keV mass is known to decay radiatively into a photon and an active 
neutrino, producing a narrow X-ray line with an energy Ec = [19]. The width of 

this decay line increases as the fifth power of sterile neutrino mass and as square of its 
mixing angle, being potentially detectable in various astrophysical DM sources like X-ray 
background, galaxy clusters, nearby galaxies and our own Milky Way [25]. Bounds on sterile 
neutrino radiative decay have been used to place upper limits on sterile neutrino mass. Strong 
upper limit constraints in the NRP scenario are obtained from the observations of Andromeda 
galaxy (M31) by XMM-Newton < 3.5 keV) [33, 34] and Chandra < 2.2 keV) 

telescopes [35]. When combined, the upper and lower limits appears to rule out the sterile 
neutrino NRP scenario. 

The DM sterile neutrino resonant production (RP) via resonant Mikheyev-Smirnov- 
Wolfenstein (MSW) [36, 37] conversion of active neutrinos to sterile neutrinos through the 
Shi-Fuller mechanism [38] in the presence of leptonic asymmetry has also been proposed 
[19, 20, 39]. Although the SM of particle physics predicts the leptonic asymmetry of the 
same order as the baryonic asymmetry, L ~ R ~ 10“^*^, there are scenarios in which much 
larger leptonic asymmetry can be produced [40-42]. 

The Tremaine-Gunn bound on sterile neutrino mass for the RP case leads to > 1 keV 
[34]. Also, the combined analysis of WMAP5 and Ly-a data indicates a bound of > 2 
keV in the frame of mixed Cold plus Warm Dark Matter (CWDM) model [43]. 

Sterile neutrino produced via RP mechanism was also forecast to have a radiative decay mode 
that could be detectable by X-ray telescopes as an emission line in spectra of X-ray clusters 
and galaxies [19, 44]. 

Other proposed mechanisms are sterile neutrino production via interactions with the inflaton 
field [45-47] and by the decay of heavy scalar singlets [48, 49]. 

Recently, few independent detections of a weak X-ray emission line at an energy of ~ 3.5 
keV have been found in a stacked XMM-Newton spectrum of 73 galaxy clusters with redshifts 
in the range 0.01 - 0.35, in Chandra ACIS-I and ACIS-S spectra of Perseus [50] and in the 
XMM-Newton spectra of Andromeda galaxy and Perseus galaxy cluster [51]. Evidences of 
this emission have been also searched in Milky Way with data from Chandra [52] and XMM- 
Newton [53, 54], in deep Suzaku X-ray spectra of the central regions of Perseus, Coma, Virgo 
and Ophiuchus clusters [55, 56], in the stacked XMM-Newton spectra of dwarf spheroidal 
galaxies in the vicinity of Milky Way [57] and in the stacked spectra of a sample of galaxies 
selected from Chandra and XMM-Newton public archives[58]. Until now, none of these 
searches were able to establish precisely if the origin of this line is is a signature of decaying 
dark matter [51] or it is related to any known atomic transition in thermal plasma, with 
special interest in Potassium and Chlorine atomic transitions [52-54]. The evidence (or lack 
of evidence in some cases) of this emission line in different astrophysical sites placed strong 
constraints on both hypothesis, raising several controversies on this subject [59-61]. 
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Sterile neutrinos produced through a resonant mechanism remains however the simplest 
model for DM origin of the undefined ~3.5 keV of X-ray line. The RP sterile neutrino 
parameters required to produce this signal have been recently inferred from the linear large 
scale structure constraints to produce full DM density [62] . These parameters are consistent 
with the Local Group and high-z galaxy count constraints, fulfilling previously determined 
requirements to successfully solve “missing satellite problem” and “too-big-to-fail problem”. 

The cosmological constraints on subdominant RP sterile neutrino parameters have not 
yet derived. The goal of this paper is to place constraints on RP sterile neutrino parameters 
in a ACWDM model containing a mixture of WDM in the form of RP sterile neutrino and 
CDM, by using most of the present cosmological measurements. 

In particular, the reconstructed gravitational leasing potential power spectrum obtained by 
the Planck satellite has impact on cosmological parameter degeneracies when DM sterile 
neutrino scenario is considered [63], offering several advantages against other methods [2, 64]. 
As the gravitational leasing effect depends on the DM distribution in the Universe, no as¬ 
sumption on light-to-mass bias is required. The projected gravitational leasing potential is 
sensitive to the matter distribution out to high redshifts, preventing from non-linear correc¬ 
tions required only at very small scales. In addition, unlike the galaxy clustering and the 
Ly-a forest, the projected gravitational leasing potential probes a larger range of angular 
scales, most of the signature coming from large scales. Moreover, the preference of present 
cosmological data for a small but significant leptonic asymmetry have been recently found 
[65]. 

This paper is organized as follows. In Sec. 2 we briefly review the DM sterile neutrino 
RP calculations. In Sec. 3 we describe the model, methods used in our analysis and the 
datasets we consider. We present our results in Sec. 4 where we examine the consistency 
and cosmological implications of DM sterile neutrino RP scenario. In Sec. 5 we draw our 
conclusions. 

2 Sterile neutrino resonant production calculations 

The DM sterile neutrinos with rest mass mg > 1 keV are produced through neutrino oscil¬ 
lations at temperatures close to the QCD phase transition [19, 20]. The oscillation process 
take place because the neutrino mass eigenstate components propagate differently as they 
have different energies, momenta and masses. Eigenstates of neutrino interaction include the 
active neutrinos/antineutrinos, Vajva (a=e, /i, r), which are created and destroyed in the 
SM of particle physics by weak interactions as well as the sterile neutrinos/antineutrinos, 
Vsjvg^ which do not participate in weak interactions. 

For the purpose of this work we consider three active neutrino flavors Val^a and one sterile 
neutrino with the production channel Vgjvs- In the early Universe, this con¬ 

version can occur through the medium-enhanced MSW resonant mechanism in presence of 
an initial neutrino lepton number residing in the active neutrinos [19, 20, 38]: 

r - np^) 

where riu^/np^ and are the active neutrino/antineutrino and photon number densities 
respectively, = 2C{3)T ^and (/(3) is Riemann zeta function of 3. 

Above the neutrino weak decoupling ~ 1 MeV), the weak interaction rates of neutrino 
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scattering are much faster than the Hubble expansion rate and the neutrino energy distri¬ 
bution can reshuffle into a thermal Fermi-Dirac spectrum. After T^ec, when active neutrinos 
decouple from the thermal plasma, the conversion process can be resonantly enhanced [ 66 ], 
affecting the primordial nucleosynthesis (T ~ 0.1 MeV) and the photon decoupling (T ~ 0.2 
eV) epochs [67, 68 ]. 

In this section we follow the detailed resonant production calculations first performed in Refs. 
[19, 66 ] for adiabatic and continuously sweeping active-sterile resonance, extended then to 
the non-adiabatically sweeping case in Refs. [69, 70] for the computation of the primordial 
nucleosynthesis abundance yields (the full solution) and examine how the sterile neutrino 
resonant production affects the photon decoupling. 


T 


-4 
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2.1 Lepton number depletion and evolution of resonance energies 

The MSW condition for the resonant scaled neutrino momentum Cres = p/T\res is given by: 

f 6 m? cos 26 

UV2C(3)7r-2G^£, 

where ~ is the difference of the squares of sterile neutrino and active 

neutrino eigenvalues, T is the photon/plasma temperature, Gp is the Fermi constant, 9 is 
the vacuum mixing angle and is the potential lepton number corresponding to an active 
neutrino flavor a: 

, /3 = (e, /i, r). (2-3) 

The leptonic asymmetry is most conveniently measured by the neutrino degeneracy parameter 
[71-73] defined as = Pua/Tu, where is the neutrino chemical potential and Ti, is the 
present temperature of the neutrino background [Tiy/T-^ = (4/11)^/^]. The relation between 
the lepton number of neutrino flavor a and the corresponding chemical potential is given by: 

The evolution of the neutrino flavor asymmetries with the temperature is driven by the active 
neutrino oscillation mixing parameters [74, 75]. These parameters are accurately measured 
with the exception of one mixing angle, @ 13 , that only recently started to be constrained. 
Actually, the value of 613 plays a decisive role in determining the temperature at which 
the equilibrium among active neutrino flavor asymmetries is reached [76, 77]. For values of 
013 closed to the present experimental upper bound, it is shown that the active neutrino 
oscillations are effective at a temperature T ^ 8 MeV [78] and therefore the active neutrinos 
reach the chemical equilibrium before BBN. For this reason we consider in the present work 
the same degeneracy parameter for all neutrino flavors: = ^u-r- 

As the universe expands, the temperature falls causing the resonance to sweep from low to 
higher values of the resonant scaled neutrino momentum eres- This process converts VajPa 
to Us-, reducing the potential lepton number As discussed in Refs. [19, 38], the evolution 
of C,y^ is dictated by the resonance sweep rate and the dimensionless adiabaticity parameter 
7 given by: 

-1 


7 


{6w?Yl^mpiG^P^ sin^ 20 
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where nipi is the Planck mass, H is the Hubble expansion rate and g is the total statisti¬ 
cal weight for relativistic species in the early Universe. Large values of 7 (7 > 1) result 
when the oscillation length is smaller than the resonance width, leading to efficient flavor 
transformation. If the time rate change of is larger than the Hubble expansion rate, 
the evolution is non-adiabatic with 7 < 1. The likelihood probability of a neutrino at reso¬ 
nance to make the jump between the mass eigenstates is given by Landau-Zerner probability, 
Plz = exp(— 7 r 7 / 2 ) [79]. The evolution of the potential lepton number when the resonant 
active neutrino/antineutrino momentum sweeps from 0 to tres is then given by: 

1 / 7 ^ \ ^ r^res 

= (l-e-'-'7(/-.+AJ<ix, (2.6) 


where /jy^ and are the initial neutrino/antineutrino Fermi-Dirac distribution functions: 


fuM) = 


F 2 {U) + 1 ’ 


fuc (e) 




(2.7) 


and F2{^^ua) the relativistic Fermi integrals of second order: 

F2{±u)^ ( 2 . 8 ) 

Jq + 1 

2.2 Phase-space distributions 

The time evolution of the sterile neutrino phase-space distribution, , can be described by 
the classical Boltzmann equation [19, 20]: 

^/i 2 .(e) -H{t)e ~ r(i/„ ^ J^s)[fuAe) - /^,(e)], (2.9) 

where: fiy^ is the active neutrino phase-space distribution, H(t) is the Hubble expansion rate 
and r(r'a —)■ Vs) is the sterile neutrino effective production rate: 


T{Ua -t 12s) ^ ^T^^{p,T) {Pm{l2a ^s)) , 

where Tjy^ is the collision rate: 


( 2 . 10 ) 


r..(p,r) 


1.27G‘^ppT^, a = e 
0.92G^ppT^, a = p,T 


The average oscillation probability {Pm{i2a —t i^s)) is given by: 


1 A^(p)sin^20 

2 /SP‘{p) sin^ 26 + D‘^{p) -|- [A(p) cos 26 + 


( 2 . 11 ) 


where A(p) = 6 m‘^(2p is the vacuum oscillation factor, D{p) = riy^(p)/2 is the quantum 
damping rate, V'^ is the thermal potential and is the asymmetric lepton potential: 


= 2V2C{3)TT-'^GFT^Ca ■ 


( 2 . 12 ) 


For temperatures characteristic to the post weak decoupling era (T < 3 MeV) the contribution 
of the thermal potential V'^ in Eq. (2.11) is very small and can be neglected. 
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Starting from thermal distributions at a scale factor o/oo rv-' 10 (oo =1 today), 

corresponding to a temperature T ~ 133MeV(ms/lkeV)^/^ (near sterile neutrino maximum 
rate production temperature), we simultaneously evolve Eqs. (2.2), (2.5), (2.6) and (2.9) to 
obtain the evolution of Ci,^{eres) in the expanding Universe for the entire range of resonant 
scaled neutrino momentum values. In this calculation we neglect the effects of the QCD 
quark-gluon phase transition that may disturb the DM sterile neutrino production rate by 
the alteration of the time-temperature relationship in the early Universe (due to variation of 
the statistical weight in different particles) [19, 63]. 

Figure 1 presents the dependence of tresC^vX^res) on sterile neutrino masses and mixing 
angles for two different initial lepton numbers, L 4 =12 and 24 (we define L 4 = lO^L^X- The 
plots presented show that the maximum on the ^res^vX^res) curves occurs at different scaled 
neutrino momentum values emax- For each case, the resonance sweeps adiabatically and 
continuously up to emax ■ For eres > ^max the resonance sweeps non-adiabatically resulting in 
an incomplete conversion Va and a partial depletion of . 

Figure 2 presents the dependence of the final phase-space distribution of sterile neutrino 
fusX) the comoving momentum e = p/T obtained for Vf, —)• flavor transformation. The 
models with fixed rus and sin^ 29 and increasing L 4 values have increasing average of the 
comoving momenta (p/T), leading to larger-scale cutoffs in the gravitational potential and 
matter power spectra. The same behavior is present for models with fixed L 4 and sin^ 29 and 
increasing nris values. The models with fixed L 4 and mg and increasing sin^ 29 values have 
opposite behaviors. Larger values of sin^ 29 lead to larger sterile neutrino production rates. 
The resonance through momentum distribution occurs in this case at higher temperatures 
and smaller comoving momenta, as can be seen from Fqs.(2.2). 

For this computation we assume a background cosmology consistent with the most recent 
cosmological measurements [2] with energy densities of Qm = 0.28 in matter, Da = 0.72 
in cosmological constant, a Hubble parameter of Hq= 69 km s~^Mpc~^ and a total active 
neutrino mass m,y=0.35 eV, as indicated by the minimal extension of the base ACDM model 
by the addition of three degenerated active neutrino flavors (see model Ml below). 

2.3 Effects of DM sterile neutrino resonant production on ^eff 

We examine the impact of DM sterile neutrino resonant production on the effective number 
of relativistic degrees of freedom, A^e//) parameter related to the radiation energy density at 
the photon decoupling. We infer Neff from the value of r^/r^, the ratio of the sound horizon 
rg to the photon diffusion length at the photon decoupling, parameter proposed to be used 
as a proxy for W// in the case of cosmological models that involve non-equilibrium neutrino 
distributions [67, 68 ]. 

A simple scaling analysis [67] show that a larger Hubble rate results in a smaller value of the 
ratio rg/rd and since A{rg/rd) ~ —ANgff {vg/rd is monotonically decreasing with Ngff) this 
imply that a larger Hubble rate leads to a larger value of Neff- 

Employing the same background cosmology as in previous subsection and a primordial helium 
mass fraction Yp=0.24, we compute the evolution with the scale factor of rg and given by 
the following integrals: 


a'^77(a') 3(1-|-7?) ’ Jo a'‘^H{a') ane{a')aT 6(1-|-7?)^ 
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Figure 1. Left: The dependence of eres^i/ei^res) on sterile neutrino mass TOs= 7 keV (continuous 
line) and ms=2 keV (dashed line) for initial lepton numbers L 4 = 12, 24 and sterile neutrino mixing 
angle sin^ 20 = 7 x 10”^^. Right: The dependence of tres^^vA^res) on sterile neutrino mixing angle 
sin^ 20=7 xl0“^^ (continuous line) and sin^ 20 = 7 x 10“^° (dashed line) for initial lepton numbers 
T 4 = 12,24 and sterile neutrino mass ms=7 keV. We define L 4 = lO^Lj^^. 


where H{a') is the Hubble expansion rate, nAa') is the free-electron number density, ax is 
the Thomson cross section and R{a) = 3 /? 6 /( 4 / 9 ^) is the ratio of the baryon energy density 
yOfo, to the photon energy density p^. We also use the recombination code Recfast [80] to 
evolve the free electron fraction Xg and the number density of free electrons Ue with the scale 
factor, required for the computation of 

In Figure 3 we plot the evolution with the scale factor of the changes in N^ff and Xg relative 
to the model with zero lepton number for sterile neutrino resonant production models with 
L 4 = 4.3 and 7.5 (RP-case) and for the leptonic asymmetric models with the same values of 
L 4 (Li/-only case). We see a larger value of Xg// at the photon deconpling (that imply a 
higher Hubble expansion rate) in both cases, relative to the case with zero lepton number. 
The increase of Xg// is smaller in the RP case relative to the Ly-orAy case with the same of 
lepton number. 

Also, a larger Hubble expansion rate at the photon decoupling implies a larger value of the 
free-electron fraction Xg and an earlier epoch of Xg freeze-out. 

3 Cosmological constraints 
3.1 Method and datasets 

We compute the angular power spectra of the CMB anisotropies, the matter density flnc- 
tuations power spectrum and the gravitational lensing potential power spectrum with the 
Boltzmann Code for Anisotropies in the Microwave Background, CAMB [81], modified to 
allow the resonant flavor conversion production calculations as presented in the previous 
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Figure 2. The dependence of the final phase-space distribution of sterile neutrino on the comoving 
momentum e = p/T obtained for —>■ Vs flavor transformation. From top to bootom the plots 
represent sterile neutrino phase-space distributions of: ms=7 keV and sin^ 20 = 7 x 10“^^ models 
with L4=2.5,4.3,7.5 and 8.7 (left panel), L4= 4.3 and sin^ 20 = 7 x 10“^^ models with ms=4,7,10 and 
15 keV (middle panel) and ^4= 4.3 and ms=7 keV models with sin^ 20= 3 xl0“^°, 1.5 xl0“^° and 
7 X 10”^^. For each case we also indicate the sterile neutrino energy density parameter, flgh^ at the 
present time. 


section. The Vajva and Vg phase-space distributions alter the homogeneous Vajva and en¬ 
ergy density and pressure, the perturbed inhomogeneous components and the gravitational 
source term in the Boltzmann equation that depends on the logarithmic derivative of the 
phase-space distributions with respect to comoving momentum. 

The impact of the leptonic asymmetry on cosmological observables has been discussed in the 
literature [82-84]. In particular, the non-zero leptonic asymmetry change the neutrino veloc¬ 
ity dispersion, the neutrino free-streaming length and the Jeans mass [85, 86]. The leptonic 
asymmetry also increases the radiation energy density, usually parametrized by variation of 
the number of relativistic degrees of freedom AA^e//(Ci^) given by: 


= 3 




(3.1) 


Figure 4 presents the gravitational lensing potential power spectrum and the matter 
density fluctuations power spectrum T’(A:) for = 7.14 keV, sin^ 20=6.46 x 10“^^ and lepton 
asymmetry numbers L4=2.5,4.3 and 7.5, having larger cutoff scales. 

The leptonic asymmetry shifts the beta equilibrium between protons and neutrons at the BBN 
epoch, leading to indirect effects on the CMB anisotropy through the primordial helium mass 
fraction Yp that decreases monotonically with increasing ^e- 

The primordial helium mass fraction is also sensitive to the shape of phase-space 

distribution. As the phase-space distributions determine the rates of the neutron and 
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Figure 3. Evolution with the scale factor of the changes in the effective number of relativistic 
degrees of freedom, N^ff, (left) in and of the free electron fraction, (right) for two lepton number 
values, ^4=4.3 (continuous lines) and ^4=7.5 (dashed lines), relative to the case with zero lepton 
number, for sterile neutrino resonant production models (red lines - RP) and and leptonic asymmetric 
models (blue lines - only). The vertical dotted line corresponds to the epoch of photon decoupling 
(adc 10 - 4 ). 

proton interaction rates at BBN, the t'e/^e spectra will change these rates and hence the 
helium mass fraction over the case with thermal Fermi-Dirac spectrum [66, 69, 70]. We use 
the BBN PArthENoPE code [87, 88] to compute the helium abundance, Yp, for different 
values of AN^ff, and changes of neutron and proton interaction rates. 

We split the total variation of the radiation energy density in two uncorrelated contributions: 

AN,fj = AN,jf{i,) + AN°ff , (3.2) 

first due to the net leptonic asymmetry of the neutrino background as given by Eq. (3.1) 
and second due to possible extra contributions from unknown processes. 

We adapt the latest version of the publicly available package CosmoMC [89] for our 
cosmological analysis and use the following datasets and likelihood codes: 

• The Planck CMB anisotropy angular power spectrum, combined with WMAP-9 
year polarization power spectrum at low t [90] and the corresponding codes [1, 91]: 
Commander, that computes the low-/ Planck likelihood, CamSpec, that computes the 
Planck likelihood for 50 < / < 2500, LowLike, that computes the likelihoods for 
2 < / < 32 polarization data and Lensing, that computes the likelihoods from Planck 
lensing power spectrum for 40 < / < 400 [92]. 

• The high-/ CMB data from Atacama Cosmology Telescope(ACT) [93, 94] and the South 
Pole Telescope (SPT) [95, 96]. 
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Figure 4. The reconstructed gravitational lensing potential power spectrum Cf’'^ (left) and the 
matter density fluctuations power spectrum P{k) (right) for mg = 7.14 keV, sin^ 20=6.46 x 10“^^ and 
lepton asymmetry numbers L4=2.5,4.3 and 7.5 having larger cutoff scales. 


• The baryon acoustic oscillation (BAO) from the Sloan Digital Sky Survey (SDSS) Data 
Release 7 (DR7) [97] at redshifts z= 0.2,0.35, the reanalyzed SDSS DR7 galaxy catalog 
data at z= 0.35 [98, 99], the SDSS Baryon Oscillation Spectroscopic Survey(BOSS) 
Data Release 9 (DR9) at z = 0.57 [100, 101] and the 6dF Galaxy Survey (OdFGS) [102] 
BAO data at z=0.1 [102, 103]. 

3.2 Analysis 

In this section we evaluate the impact of sterile neutrino resonant production mechanism on 
the cosmological parameters. We consider the following extensions of the standard ACDM 
cosmological model: 

• Ml: ACDM+Smj/^: 

The minimal extension of the base ACDM model by the addition of three degenerated 
species of massive neutrinos with total mass (within which ACDM is nested at 

Neff=2>MQ and yp=0.24). 

• M2: ACDM+Smj/^ +S^jy^+AA/e// +Fp- 

The addition to Ml of the neutrino chemical potential the extra radiation energy 
density AAg// given in Eq. (3.2) and the BBN prediction of the primordial helium 
abundance Yp. 

• 1VI3: ACDlVI-bXjTTrt/^ “bE^^^-bAAgyj TTig-\-sui^ ‘20-\-Yp'. 

The addition to M2 of one DM sterile neutrino with the mass m*, the matter mixing 
angle sin^ 29 and the BBN prediction of the primordial helium abundance Yp. 

Our main results are summarized in Table 1. 
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Table 1. The table shows the mean values and the absolute errors on the main cosmological param¬ 
eters obtained from the fits of different extensions of the ACDM model discussed in the text with 
PLANCK-|-WP-|-highL-fBAO-flensing dataset. For all parameters, except and sin^ 26 we quote 
the errors at 68% CL. For Em^ and sin^ 26 we give the values of 95% upper limits. 


Parameter 

Ml 

M2 

M3 

Qbh‘^ 

0.02227±0.00019 

0.02221±0.00024 

0.02209±0.00023 

^cdm^ 

0.1161±0.0015 

0.1162±0.0015 

0.1169 ±0.0021 

1OO0MC 

1.04156±0.00059 

1.04222±0.00095 

1.04261±0.00072 

r 

0.0882±0.0091 

0.0851±0.0092 

0.0814±0.0121 

^rriu^ (eV) 

< 0.35 

< 0.25 

< 0.21 

ms(keV) 

- 

- 

6.08 ±1.62 

sin^ 26 

- 

- 

< 5.61 xl0-i° 

^Va 

- 

-0.062±0.079 

0.000002±0.000093 

^eff 

- 

3.37±0.19 

3.25±0.15 

rig 

0.9586±0.0046 

0.9679±0.0048 

0.9711±0.0052 

lnl0i°%^ 

3.213±0.019 

3.231±0.023 

3.192±0.041 

Da 

0.699±0.012 

0.712±0.011 

0.713±0.012 


0.302±0.012 

0.288±0.011 

0.287±0.011 

0-8 

0.793±0.021 

0.801±0.022 

0.772±0.015 

Ho 

68.16±1.02 

70.15±1.36 

69.95±1.17 

Yp 

0.24 

0.267±0.019 

0.24968±0.00011 

Zeg 

3306±33 

3221 ±36 

3183±46 

Lua 

- 

-0.015±0.012 

0.000123±0.000022 

rsird 

6.478±0.011 

6.389±0.032 

7.013±0.011 


3.2.1 Neutrino masses and mixing 

Top left panel from Figure 5 presents the likelihood posterior distributions of total active 
neutrino mass, showing a decrease of the upper limit (at 95% CL), from 0.35 eV, for 

the case of minimal extension of the ACDM model (Ml), to 0.25 eV when only the leptonic 
asymmetry is included (M2) and a subsequent decrease to 0.21 eV for the RP sterile neutrino 
scenario (M3). 

In the top right panel of Figure 5 we plot the likelihood posterior probability distribution of 
the sterile neutrino mass fraction = U.gh?■ The dominant effect on sterile neutrino 
resonant production is given by the value of the initial potential lepton number, rather 

than sterile neutrino mass. This is shown in the bottom panels from Figure 5 that present 
the joint conhdence regions £(((** - and mg - fu^ ■ We find < 0.078 for a value of lepton 
number per flavor L4 = 1.23 ± 0.04, sterile neutrino mass rus =6.08 ± 3.22 keV and mixing 
angle sin^ 20 <5.61 xlO”^*^ (errors at 95% CL). 

Figure 6 presents the marginalized likelihood posterior distribution obtained for and the 
joint confidence regions mg -sin^ 26. 

3.2.2 ^He abundance yield 

The BBN '^He yield depends sensitively t'e/T'e degeneracy parameter, as well as on the 
shape of its phase-space distribution function. We find the preference of cosmological data 
for smaller values of and ^He mass fraction, Yp, in the case of RP sterile neutrino scenario 
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Figure 5. Top: The marginalized likelihood posterior distributions for active neutrinos total mass 
and and sterile neutrino mass fraction, as obtained from the fits of different 

extensions of the ACDM cosmological model with PLANCK+WP+highL+BAO+lensing dataset (see 
also the text). Bottom: The joint confidence regions /I™* - and rris - (at 68% CL and 95% 
CL). 


(M3) when compared to the similar predictions obtained in the model with leptonic asymme¬ 
try only (M2). This is shown in Figure 7 where we plot the marginalized likelihood probability 
distribution of obtained from the ht of model M3 and the joint conhdence regions - 
Yp from the fit of model M2 and M3 with our dataset. We obtain ^i,=0.0005±0.0004 and 
yp=0.2497 ± 0.0002 (95% CL) for M3 case, in very good agreement with the Planck deter¬ 
mination [2], Yp = 0.2485±0.0002, and the latest observational bounds yp=0.24968±0.00022 
[104]. 

4 Conclusions 

Recently, few independent detections of a weak X-ray emission line at an energy of ~ 3.5 keV 
seen toward a number of astrophysical sites have been reported. Until now, none of these 
searches were able to establish precisely the origin of this line. If confirmed, this signal could 
be the signature of decaying DM sterile neutrino with a mass of ~ 7.1 keV. 

The DM sterile neutrino production via resonant MSW conversion of active neutrinos to 
sterile neutrinos through the Shi-Fuller mechanism [38] is the simplest model for the DM 
origin of the undehned ~3.5 keV of X-ray line. 
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Figure 6. Top: The marginalized likelihood posterior distributions obtained for the sterile neutrino 
mass, rris and the mixing angle sin^ 26. Bottom: The joint confidence regions -sin^ 26 (at 68% CL 
and 95% CL). 




Figure 7. Left: The marginalized likelihood probability distribution of from the fit of model M3 
with PLANCK+WP+highL+BAO+lensing dataset. Right: The joint confidence regions - Yp (at 
68% CL and 95% CL) from the ht of model M2 (black) and M3 (red) with the same dataset. The 
latest observational bounds [104] are also plotted (at 68%CL). 
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Figure 8. Two-dimensional marginalized joint probability distributions at 68% and 95% CL (red 
lines) showing the dependence of RP sterile neutrino mass, ms, and the matter mixing angle, sin^ 26, 
obtained from our cosmological analysis. We also show the 95% CL upper limits reproduced from 
Ref. [32] based on deep Chandra [32] and XMM-Newton [35] observations of M31 and Suzuku observa¬ 
tions of Ursa Minor [55, 56]. The vertical lines represents the lower mass limits from Tremaine-Gunn 
phase-space considerations (m-s ~ 0.4) keV [27], Coma Berenices phase-space ~ 1.5) keV, Segue 

I phase-space ~ 2.5) keV and M31 subhalo counts ~ 8.8) keV. The upper and lower 

bounds (95% CL) of sterile neutrino mass obtained in Dodelson-Widrow model [23] are also indicated 
(green dashed lines). The yellow symbol corresponds to the detection from Ref. [50] (mg ~ 7.1 keV, 
sin^ 26 = 7 X 10“^^ and 90% statistical error bar). 


The RP sterile neutrino parameters required to produce this signal have been recently in¬ 
ferred from the linear large scale structure constraints to produce full DM density [62] . These 
parameters are consistent with the Local Group and high-z galaxy count constraints, fulfill¬ 
ing previously determined requirements to successfully solve “missing satellite problem” and 
“too-big-to-fail problem”. The cosmological constraints on subdominant RP sterile neutrino 
parameters have not yet derived. 

In this paper we place constraints on RP sterile neutrino parameters in a ACWDM 
model containing a mixture of WDM in the form of RP sterile neutrino and CDM, by using 
most of the present cosmological measurements. 

We make a coupled treatment of the weak decoupling, primordial nucleosynthesis and photon 
decoupling epochs in the sterile neutrino RP scenario including the extra radiation energy 
density parametrized by the effective number of relativistic degrees of freedom A^e//- 
We compute the radiation and matter perturbations including the full resonance sweep so¬ 
lution for ValvoL I's flavor conversion in the expanding Universe and provide constraints 
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on the cosmological parameters and sterile neutrino properties, by using most of the present 
cosmological measurements. In particular, the high precision CMB temperature anisotropy 
and the reconstructed CMB gravitational lensing potential power spectra obtained by the 
Planck satellite have impact on cosmological parameter degeneracies when DM sterile neu¬ 
trino scenario is considered. 

We find that the values of the main cosmological parameters are in agreement with the pre¬ 
dictions of the minimal extension of the base ACDM model except for the active neutrino 
total mass that is decreased to <0.21 eV (at 95% CL). 

Figure 8 presents the - sin^ 26 joint probability distributions obtained from our cosmo¬ 
logical analysis, compared with the most of the existing similar constrains. The cosmological 
measurements are in agreement with the sterile neutrino RP scenario with the following pa¬ 
rameters (errors at 95%CL): rus =6.28 ± 3.2 keV, mixing angle sin^ 29 <5.61 xl0“^^, lepton 
number per flavor L 4 = 1.23 ± 0.04 and sterile neutrino mass fraction < 0.078. 

Our results are in agreement with the sterile neutrino resonant production parameters in¬ 
ferred in Ref. [62] from the linear large scale structure constraints to produce full Dark 
Matter density. 

This reflects the sensitivity of the high precision CMB observables to the physics that de¬ 
termines the neutrino energy spectra at different epochs, that require a self-consistent and 
coupled treatment [68]. 
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